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Abstract 

The linear sigma model with pions and sigma coupled to baryons and a (massive) gauge vector 
meson is considered for the description of a finite baryonic density system. The stability equation 
for the bound system is studied in the homogeneous limit considering quantum fluctuations for the 
scalar and pseudo scalar fields within a variational method with truncations. Nearly exact solutions for 
the stability equation which also solve the field equations are sought semi-analytically. Bosonic fields 
are found to have non zero expected classical values at finite density corresponding to condensates of 
dynamical symmetry breakings. The two components of baryon field, neutron and proton for example, 
acquire different masses due to the isospin symmetry breaking and these states may even oscillate in the 
medium. There is a non trivial solution for the (massive) vector meson field indicating the occurrence 
of a gauge dynamical symmetry breaking at finite density typical of a superconductor. 
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1 Introduction 



When a theory, given by a Lagrangian, is invariant under a transformation group and the ground state 
of this theory is not, this symmetry is spontaneously (or dynamically) broken (SSB) [j]. 

Quantum Chromodynamics (QCD) is believed to be the fundamental theory for the Strong interac- 
tions, at least at a certain level. It has a very complex non abelian structure and strong coupling constants 
for processes at low energies. At not high energies quarks and gluons, whose charge is the SU(3) color, 
are confined and no colored state can be observed. Due to these features it is very difficult to obtain exact 
solutions. Descriptions of low energies processes and systems are mostly accomplished in finite lattices 
(at zero densities) where space-time is discretized. One is therefore lead to construct effective models 
which respect the main properties and symmetries of the QCD for the range of energy and densities for 
the hadronic process of interest. It is interesting to note that one may have different effective theories 
for certain ranges of these variables which can be different if one consider different systems such as the 
vacuum or high energy density matter such as formed in relativistic heavy ions collisions in RHIC and 
CERN. In the vacuum, the lightest strong interacting particles are known to respect, at least approxi- 
matedly, chiral symmetry whose group is the SUl(2) x SUr(2) which is spontaneously broken down to 
SU(2). This means that the QCD vacuum has a preferencial direction in chiral space corresponding to 
a scalar state. Pseudo scalar pions have small masses in the hadronic scale and they are viewed as the 
Goldstone bosons of such SSB. The vacuum acquires a non trivial structure due to the formation of scalar 
quark - anti-quark condensate < qq >, the order parameter of the Chiral SSB. Due to this condensate, 
vacuum is not invariant under chiral transformations. These features can be taken into account via sigma 
models which, in the linear realization [^, pj with mesons, implement chiral symmetry with two fields: 
the (pseudo-scalars) pions and the (scalar) sigma - which is called the chiral partner of the pion. The 
expected value of sigma is, in this model, the order parameter of chiral SSB, < a >oc< qq >. Long time 
ago S. Weinberg pointed difficulties of the linear sigma model, which are being rediscussed, for describing 
hadronic processes and proposed the non linear realization of chiral symmetry in which the sigma field is 
"eliminated" ||. There are nowadays strong evidences for the existence of the sigma || although there 
is no consent j|] . Both implementations of chiral symmetry have been extensively developped since then. 

The linear realization of chiral symmetry exhibits advantadges over the non linear realization mainly 
for the description of dense systems ||. Pions are expected to have zero classical components (conden- 
sates) < 7r >vac= in the vacuum. At finite density, QCD is known and expected to have a very complex 
phase diagram with the appearance of other condensates |M] with the possibility of a restoration of the 
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chiral symmetry when the chiral condensate < qq > is expected to be close to zero at high energy den- 
sities. With different approaches, finite density QCD and effective models have been intensively studied 
in the last years. The behavior of the symmetry order parameter (s) with density (and temperature) is a 
very relevant aspect for the understanding of matter and its fundamental interactions as well as for the 
description of the evolution of the Universe. 

In this work we study the chiral Linear Sigma Model (LSM) with baryons and a (massive) vector 
meson. The mass of the vector meson is considered to be non zero only when the baryonic density is 
different from zero. The exact field equations are truncated to allow for analytical solutions by considering 
particular prescriptions for the stability condition of the system. We found that there is a (isospin) 
spontaneously broken symmetry generating non zero expected value for the pion field at finite densities 
(not necessarily high) . The quantum fluctuations of the sigma and pion fields are considered in a truncated 
version of the variational approach with Gaussian wave functionals. The solutions for the bosonic fields 
are discussed and some other possible consequences are analyzed. 

2 Linear sigma model at finite density 

The Lagrangian density of the Linear Sigma Model for nucleons ./V(x), sigma and pions (a, tt) covariantly 
coupled to a massive gauge vector meson is given by: 



where the covariant derivative is: = <9 M — igyV^, the gauge invariant tensor is: = d^V" — d v V iX . 
9Vi 9s and A are the coupling constants and v = /„■, the pion decay constant, in the vacuum. The 
gauge invariant coupling of the temporal component of the vector meson to the baryons is equivalent 
to a redefinition of a chemical potential. Pion mass is known to break the chiral symmetry explicitly, 
in such case one should add another (small) term C s h = ca, where c is a constant proportional to the 
mass. This would be of no relevance for our main results. We consider the possible existence of classical 
expected values (condensates) for all the mesonic fields and look for solutions for the equations of motion. 
In this work the whole baryon masses come from the the coupling to the scalar mesonic field by a Higgs 
mechanism (M* = gs&)- 

The vector meson mass rriy is considered to be zero at zero baryonic density (ps) and non zero at 
Pb / 0. We will be concerned, in this level of analysis, with an homogenous system in which case the 
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other (spatial) components of the vector meson Vi would in part produce small shifts in the numerical 
results. It would only appear in the vector meson mass term at finite density and in the term of coupling 
to baryons. A further account of a complete coupling will be given elsewhere. 

The baryon (nucleon) field is quantized in terms of creation and annihilation operators. Its wave func- 
tion can be written as superposition of spinor (x), isospinor (77) and coordinate components u(p), v(p). 
It generates non zero scalar, baryonic and pseudo-scalar densities (p s , pp and pps)- We will not explicitly 
evaluate here all these quantities but only make use of pb and the quantized fermionic energy density 
pj, given respectively by: 

In these expressions kp is the baryon momentum at the Fermi surface. 

To take into account the quantum fluctuations of the sigma and pion fields we consider a truncated 
version of the variational approach using a Gaussian trial wave- functional in the Schroedinger picture || . 
With the variational principle an estimation for the energy density of the ground state can be obtained by 
calculating the averaged energy with trial wave- functions (Tl^ =< &\H\& >) whose (trial) parameters are 
fixed when the energy density is minimum with relation to its parameters. Although this statement is not 
very precise due to the subtraction procedure required by renormalization for the ultraviolet divergences 
one can not neglect such idea. The averaged value of the Hamiltonian is calculated with trial Gaussian 
wave-functionals for the scalar and pseudo-scalar fields: < ^[cr, tt]\H(ct, 7?) | ^ [ct, tt] >. For the sigma we 
can write: 

* [<r(x)] = Nexp^~ J dxdy^(x)G s 1 (x,y)^(y)| , (3) 

Where 5a(x) = c(x) — a; the normalization is N, the variational parameters are the condensate 
a =< ^'|cr|^' >, the quantum fluctuations represented by the width of the Gaussian Gs(x,y) =< 
^\a(x)a(y)\^/ >. An analogous expression for the pion sector is considered with variational parameters 
given by: tt =< tt > and G p ,b , which is a matrix that can be considered to be diagonal as a particular 
case along this work {G a p a = Gp). This reduces the corresponding functional space but it guarantees the 
explicit chiral invariance. We will assume that these quantum fluctuations only intervene for the meson 
masses as well as for a shift of the respective fields: 

~i = 7? 2 + G P , a 2 = a 2 + G s . (4) 

This corresponds to a truncation on the self consistency with a particular renormalization energy point. 
This approximation corresponds to neglect terms in the energy density calculated with the Gaussian 
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functionals. If all terms are considered we would only have redefinition of the energy density. It does not 
need self consistency to be fixed. More exact calculations are being done for this model in ||. 

The minimizations of the averaged energy with respect to the Gaussian variational parameters yield 
the GAP and condensate equations which define the minimum of the potential for these fields. The 
following set of equations is obtained for the sigma sector: 

6 -^ = 0^a\(a 2 + 3G s + t 2 + G P -v 2 )+^ = 0; 

° a -2 da (5) 

/ ~ 2 

From this last expression we can write the following expression for the sigma mass: fj/g = A ( 3a 2 + tt — v 2 + 
For the sake of clearness we truncate the exact solution considering that: 1) fermionic density does not 
depend on the fluctuations G{, 2) fluctuations are small and equal {Gp = Gs)- The corresponding 
equations for the pion field are: 

\7t a (t 2 + 2G P + d 2 -v 2 ) + ^ = 



Hp 



d^a (6) 



A (3^ + a 2 - v 2 + ...) -> 

The first of these two expressions contains one of the most relevant results of the present work. As far as 
the fermionic density depends on a, via the in medium mass M*, and the equation for a (expression (|5[) 
does depend on n a we find that the pion condensate n a must be nonzero whenever the fermionic density 
(and therefore baryonic density) is not zero. We have found therefore that: 

dp f = d Pf da Q 
dir a da dit a 

As discussed above, the expressions for the meson masses p 2 s ~ 500MeV and //p ~ 0, and the shift of the 
fields (Q), are the only effects of the quantum fluctuations in the present work. Since we will not explicitly 
calculate the two point function Gi, there will be no concern with Ultra Violet divergences here. 

The Euler Lagrange equation for the vector meson was calculated for a gauge in which only the 
component Vq is non zero and homogeneous. If we consider the other components V{ 7^ the equations 
are slightly changed, but the conclusions remain valid with another coupled equation. The equation is 
given by: 

9v(pb + V ^] -m 2 / y = 0. (8) 
The total averaged energy density can be written as: 

U = Pf + g V V pB ~ \m 2 v V 2 + ^(a 2 + ~rr 2 - v 2 ) 2 . (9) 



In this work solutions for the above equations (|],|1J8|) are researched such that the main properties of 
a (bound) finite density system are consistently described. The stability condition of the system can be 
written as: 

< 0, (10) 



dp B 



Pb=Po P b Pb=Po 

where po is the stable density. To guarantee that this expression is satisfied we consider some (approx- 
imated) prescriptions for the dependence of the involved variables in on the baryonic density. Namely, 
we divide the expression for the energy density (||) for the condition ( |i~0| ) into three separated ones: 

d Pf = Pf 

dpB PB ' 

d{a 2 + # 2 - v 2 ) _ (a 2 + # 2 - v 2 ) r n s 

dp B PB 
dTiy Tiy 

dpB PB 

In this last expression Hy is the energy density with contributions of the vector meson. A detailed com- 
parison of this prescription with the exact result will be shown elsewhere. It is certain that the solutions 
for these (stability) equations do obeys the equations of motion of the fields, being this approximation 
therefore reliable. 

From the first of the differential equations (11) we find a solution for the dependence of pj on the 
baryonic density {pf = Pj(pb)) which is in excellent agreement with that resulting from the integration 
of expression (Q) ||. One more evidence that the prescriptions of expressions ( JXT] ) are good. 



From the second expression in (11) we find a constraint which defines a symmetry radius in the 
medium: 

(a 2 + ~T?-v 2 ) = C y fp B -. (12) 

In this expression C is a constant. Therefore in the vacuum: a 2 ~ v 2 = f 2 as discussed above. We stress 
that a = f v a in the vacuum. This means that the quantum fluctuations change the condensate value. 

The GAP and condensate equations (^J, |6]) , for a and for the 7r can be faced as differential equations 
for pf and they can be written as: 

[a + V ' \da 2 - \ d ^ 2 {L6) 

Where v 2 = v 2 — G. These equations are isomorfic and show an equal dependence of pf with each of 
these (chiral) fields. An approximated solution for the two condensate equations can be found as if pj 
were a function of these fields independently. This yields Pf = pf\^) and Pf = Pf (n) with constants 
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fixed for the boundary ps = where 7? = and a = v. Equating p^ (solution of the first equation) to 

the second equation we find the following approximated value for the pion condensate (if |7f 2 | << v 2 ): 

~_ 2 a 2 (a 2 -v 2 ) , x 

vf ~ — y — '- (14) 

~2 

In these solutions, as well as in others more exact, tt may be either positive or negative. It would 
correspond respectively to a real or complex condensate value. This solutions do reproduce values close 
to the exact solutions. Exact solution for the condensate equations will be given elsewhere. 

We can also add the two differential equations ( |l3|) which can be seen as partial differential equa- 
tions. We equate the solutions for these equations and obtain the following constraint expression for the 
condensates: 

^v^=^(a 2 -y) 2 -(# 2 -y) 2 ). (15) 
Where Cf is a constant from the partial differential equations. This also expresses the dynamical symme- 
try breaking which occurs in the medium mainly because CfC^/ps / 0. We obtained a consistent basis 
for the study of the behavior of the condensates (as order parameters) from a Quantum Field Theory 
which has a SSB with density. 

Finally, considering the equation of Vq - expression (g) - as a differential equation of the baryonic 
density pb as a function of Vq we obtain the following solution: 



_ _ . , -g V pB±\/gvPB- 2C VpBmv Ma\ 
V (PB) = 2 ' { } 

where Cy is a constant. This constant will be the only contribution of the vector meson sector to the 
energy density Tiy = Cyps- In the limit of zero density we get Vq — > 0. It is seen that the baryonic 
density generates a non zero value of Vo - which can be viewd as a condensate. This may be another 
dynamical symmetry breaking. Requiring the energy to be stable with relation to the solution Vb we 
find that: Cy = 2gyVo. Where Vq is the value of Vb at a stable density p$. In this point we have 
my = —2gy po/Vq. From this we see that the mass of the vector meson is proportional to the density, i.e., 
it is an in medium effect, as required. The massive vector meson is characteristic from a superconductive 
state which occurs at finite density. The question of (non) renormalizability (in a finite density medium) 
of the massive vector field will not be addressed here. 

2.1 Discussion 

Besides the appearance of another condensate, n, at finite baryonic density due to the isospin SSB we 
find the possibility, although not very appealing, of further chiral symmetry breaking with a > v at finite 
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baryonic density as found in [ 10 1 . The numerical studies will be shown in a forthcoming work. With the 
present results it is not yet possible to determine, unambiguously, the symmetry behavior with density 
in the present model as precluded in other works ||. 

We found that, at finite baryonic densities, there is a non zero pion classical expected value (conden- 
sate). Departing from the Lagrangian density or from the equation of motion of the baryons we propose 
the following definition for their effective masses. This would also explicit the feature that the nucleon 
densities depend on the pion condensate. From the averaged value of the energy we may consider a 
matrix for M* which depends on the isospin (and eventually spin) of the nucleons: 

K,b;s = 9s< *k,7r]|. < N a>s \(a + i~f 5 T.rf)\N b , s > .|*[<r,7r] >= g s {o + iM d {a ^ d ). (17) 

In this expression a,b,d stands for the isospin number (neutrons/protons), and M is a non diagonal 
isospin matrix. This allows for the possibility of different values of baryonic masses - i.e. a mass splitting 
- but it also means possible oscillations between the isospin states (proton and neutron for example) 
in the medium. This matrix includes a dependence on the baryon spins as well which is not discussed 
further here. However, the quantities calculated with it in the present work -mainly the densities (p/)- 
were the usual ones, i.e., with a constant diagonal effective mass due to a, i.e. M* b (n) ~ M* , which is 
the leading term. This non trivial solution corresponds to a non invariant ground state under an isospin 
transformation, although the Lagrangian is symmetric. The massive vector meson is characteristic from 
a superconductor - finite density. It may be the photon itself which would acquire mass in the presence 
of the hadronic field. In this case, arguments of the kind of vector meson dominance to explain its usual 
identification with the omega would be appealing. The velocity and other properties in the baryonic 
medium of this massive meson will be calculated elsewhere Q. 

Phase diagram of QCD would then include this pseudo scalar condensates. Consequences for de- 
scription of effects in (relativistic) heavy ion collisions (as well as for astrophysical dense objects) can 
be expected: baryon properties such as masses are modified, with possible splittings and baryon states 
oscillation. Pion multiplicities in (relativistic) h.i.c. may increase with the presence of < tt > with respect 
to predictions without it. The present model will be considered for the description of (hadronic) nuclear 
matter properties and eventually of atomic nuclei elsewhere as well as inhomogeneous situations where 
topological properties may arise [p. 
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